The entangling power of passive optical elements 
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We investigate the entangling capability of passive optical elements, both qualitatively and quantitatively. 
We present a general necessary and sufficient condition for the possibility of creating distillable entanglement in 
an arbitrary multi-mode Gaussian state with the help of passive optical elements, thereby establishing a general 
connection between squeezing and the entanglement that is attainable by non-squeezing operations. Special 
attention is devoted to general two-mode Gaussian states, for which we provide the optimal entangling proce- 
dure, present an explicit formula for the attainable degree of entanglement measured in terms of the logarithmic 
negativity, and discuss several practically important special cases. 



Entangled states of light field modes may be generated by 
transmitting two squeezed states through a beam splitter [1]. 
This is one of the experimentally accessible procedures for 
generating continuous-variable entanglement in optical sys- 
tems [2]. Moreover, it is a particular example of a situation 
where passive optical elements exhibit their entangling power 
when applied to Gaussian input states. It is well known that 
the presence of squeezing is necessary for obtaining entan- 
glement in this manner [1]. However, the degree of the at- 
tained entanglement is by no means the same for all input 
states: it depends to a large extent on the degree and direc- 
tion of squeezing of the incoming modes and on the specific 
properties of the beam splitter. This raises the question under 
what circumstances such an entangling procedure is optimal 
in the sense of generating states which have the maximal at- 
tainable amount of entanglement. And in general, by means 
of arbitrary passive optical elements, what are the require- 
ments such that entanglement can be generated between any 
bi-partite split of a system in a multi-mode Gaussian state? 

In this letter we address the question of the entangling 
power of passive optical elements acting on any number of 
modes in an arbitrary Gaussian state, qualitatively as well as 
quantitatively. Passive optical operations can be implemented 
by using beam splitters and phase shifters [3]. These are cheap 
operations and easy to implement in contrast to squeezing op- 
erations. Therefore we will consider squeezing as a potential 
resource for entanglement and ask for the requirements and 
the optimal way of entangling a squeezed state by means of 
passive operations, which we assume to be available in arbi- 
trary quantities. The main result and starting point is a nec- 
essary and sufficient condition for the possibility of creating 
distillable entanglement on general Gaussian initial states - 
pure or mixed - between any bi-partite split of an n-mode 
system with the help of passive optical elements. We then 
introduce a lower bound for the attainable degree of entangle- 
ment measured in terms of the logarithmic negativity [4] for 
n-mode systems. Moreover, we derive a general formula for 
the largest degree of entanglement of an arbitrary subsystem 
consisting of two modes. The operations that can be imple- 



mented with passive optical elements can be identified with 
the non-squeezing operations. In this sense we establish a 
quantitative connection between the degree of squeezing of 
a Gaussian state and the degree of entanglement that is attain- 
able with the application of non-squeezing operations. Of par- 
ticular interest is the case where only two modes are present. 
We will discuss this situation in more detail by explicitly con- 
structing the optimal entangling procedure and discussing sev- 
eral meaningful special cases. 




FIG. 1. n field modes in a Gaussian state p can be 
(NPPT)-entangled with respect to a partition into ha + tib — n 
modes by means of passive optical elements if and only if A1A2 < 1, 
where Ai and A2 are the two smallest eigenvalues of the covariance 
matrix associated with p. 

We start by introducing the formalism that we will use ex- 
tensively. Gaussian states are completely characterized by 
their first and second moments, where only the latter, given 
in terms of a covariance matrix V, carry information about en- 
tanglement and squeezing. For this reason we will set the first 
moments to zero, which can always be achieved by unitary op- 
erations on individual modes. The covariance matrix is then 
given by = 2(RkRi) — icT^i, k, I = 1, . . . , 2n, where the 
vector R = (Qi, . . . , Q n , Pi, . . . ,P n ) consists of the canon- 
ical coordinates for n modes and the symplectic matrix 
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governs the canonical commutation relations (CCR) 
[Rk,Ri] = iuki- A matrix represents an admissible covari- 
ance matrix if it satisfies the Heisenberg uncertainty relations 
r + ia > 0. Symplectic transformations V 1 — > S T TS pre- 
serve the CCR and therefore satisfy S T aS = a [5]. All 
symplectic transformations correspond to unitary Gaussian 
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operations [6] on the level of states, in the sense that the 
Gaussian character of arbitrary input states is preserved under 
such unitary operations. They can be decomposed [5,7] into 
active/nonlinear and passive/linear operations [8]. The latter 
can be implemented by using passive optical elements such as 
beam splitters and phase plates only [3], and are of the form 
T i — ► K T TK, 

Here U = X + iY (X, Y real) is any unitary matrix and 

si- — ( ln lln \ (3) 

relates real and complex representations by mapping cre- 
ation/annihilation operators to position/momentum operators 
via n(i2i, . . . , R 2n ) T = (oi, • ■ • ,a n ,a\,... , a) n ) T . Trans- 
formations of the type V i — ► K T TK with if as above will 
from now on be denoted as passive transformations. Any such 
if is both symplectic and orthogonal, i.e., K T K = 1, and the 
set of all symplectic transformations that can be implemented 
with passive optical elements form a group, the maximal com- 
pact subgroup K(n) of the group of symplectic transforma- 
tions 5p(2n,R) [5]. 

A Gaussian state is said to be squeezed if there exists a ba- 
sis in phase space such that at least one diagonal element of 
the covariance matrix is smaller than one. From now on we 
order the eigenvalues of T in non-increasing order, so that this 
implies that the smallest eigenvalue Ai of T is smaller than 
one [5]. Since every passive transformation K is orthogonal, 
it does not affect the squeezing of a state. 

Let us now turn to entanglement properties. A Gaussian 
state of a bi-partite system consisting of parts A and B with 
riA + ub = n modes is separable, i.e. unentangled between 
A and B, iff there exist covariance matrices Ta,^b for tia 
resp. us modes such that r > Ya ® I\b [9,10]. A necessary 
and for 1 x ns modes also sufficient condition for separability 
[9,1 1] is that the partial transpose of the state is positive semi- 
definite. This, in turn is equivalent to T > id, with the par- 
tially transposed symplectic matrix a = (1„ © E)a(l n E) 
and E = l nA ® (— t nB ) being the partial transposition oper- 
ator that reverses all momenta on one side. It has been shown 
that a Gaussian state is distillable, i.e., that its entanglement 
can be revealed using local operations and classical commu- 
nication, iff its partial transpose is non-positive [12]. 

Obviously, every entangled Gaussian state is squeezed 
since Ai > 1 would mean that T > 1 © 1 which in turn 
implies separability. Hence, a state can only be entangled by 
means of passive operations if it is squeezed initially. The fol- 
lowing proposition gives a necessary and sufficient condition 
for the possibility of transforming a general Gaussian state 
into a distillable one by means of passive transformations (see 
Fig. 1): 

Proposition 1 Let T be a covariance matrix corresponding 
to a Gaussian state of n modes. A passive transformation 



r i — ► r' = K T TK leading to an entangled state having a 
non-positive partial transpose with respect to a partition into 
riA + fiB =n modes exists iff 

AiA 2 < 1, (4) 

where Ai , A2 are the two smallest eigenvalues ofY. 

Proof: A Gaussian state of an n-mode system with covari- 
ance matrix V has a positive partial transpose iff all symplec- 
tic eigenvalues of the respective partially transposed covari- 
ance matrix T' = ( 1„ ® E)V ( 1„ ® E) are larger than or equal 
to one. The symplectic eigenvalues of f" are in turn equal to 
the square roots of the ordinary eigenvalues of — (T'd) 2 . The 
square of the smallest symplectic eigenvalue additionally min- 
imized over all passive transformations is thus given by 

i/:=inf inf (£|r 1/2 M T rMr 1/2 |0, (5) 

K ||||| = 1 

where M := KdK T and |.| denotes the standard vector 
norm. Hence, we have to show that inequality (4) is equiv- 
alent to v < 1. 

Since M is an antisymmetric orthogonal matrix, it maps 
any real unit vector onto the two-dimensional unit sphere of 
its orthogonal complement. The vector 

\^K,o)~Mr 1 / 2 \o-mO- 1/2 (6) 

therefore satisfies (£_\r^ 2 \C(K, £)) = 0. Inserting Eq. (6) in 
Eq. (5) we get 

" = " ^(flrifKf (K.oiriffK.f)) m 
>- „ £ „^„=. (£|r| « > «' |r|£ ' > ' (8) 

where the infimum in Eq. (8) is taken over all real unit vec- 
tors satisfying (£|r 1/,2 |£') = 0. This relaxes the require- 
ment that has to be of the form in Eq. (6) and there- 
fore leads to the lower bound. The minimum in Eq. (8) is 
now attained for vectors lying in the two-dimensional space 
corresponding to the two smallest eigenvalues Ai,A2 of T. 
Hence, |£) = cos0|Ai) + sin 1^2) for some (j> and oc 
\/A2sin<^|Ai) — y / Ai"cos^|A2). However, every (f> leads to 
the same value and we have 

v > AiA 2 , (9) 

showing that Ai A2 < 1 is indeed necessary for v < 1. 

In order to prove sufficiency, we have to show that there 
always exists a passive transformation K such that is of 
the form (6) and the inequalities (8,9) thus become equalities. 
Note that this is in turn equivalent to the statement that for ev- 
ery pair of orthogonal real unit vectors |£) _L 1 77) there is a pas- 
sive transformation K such that (rj\KdK T |£) = 1. We first 
show that the problem can be reduced to a two-mode prob- 
lem. Let |Ai) and | A2) be the eigenvectors associated with Ai 
and A 2 . Decomposing |Ai), | A2) into position and momentum 
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components, one may define the complex form of |Ai) and 
|A 2 ) according to 



!*<>:= ixf^+i^), , = 1,2. 



(10) 



Then there always exists a unitary U such that the vectors 
U\^>i) have only two non-zero components in the first two 
entries of the vector. In turn, according to Eqs. (2,3) this 
implies that there exists a passive transformation £ such that 
T' := S T TS has the property that the leading principle sub- 
matrix of T' associated with the first two modes has the same 
two smallest eigenvalues Ai, A 2 . 

Similarly, |£) , |?7) can be decomposed into position and mo- 
mentum components, and define |^) :— + and 
|$) := |f?W)) + t^iP)) such that = Q|f «) © ^} = 
-i=|* and analogous for \r)} and |$). Then ||*|| = 



1 and 



(r)\K5K T \g) = lm[(^\UEU^\^}}, 
<»?|0 = Re[<*|*>] = 0. 



(11) 
(12) 



Without loss of generality we fix |W) = (1, 0) T , which can al- 
ways be achieved by applying an additional unitary. Then, ev- 
ery two-dimensional unit vector |$) for which Re[($|W)] = 
is of the form 



|$) = -i(cos(2 7 ),e 2ra sin(2 7 )) 



Choosing 



U 



* cos(7) 
: sin(7) 



%a sin(7) 
: cos(7) 



(13) 



(14) 



we obtain with E = diag(l, -1) 

1 = lm[(<S>\UEU^\t>)} = (r 1 \KaK T \0, 

which completes the proof. ■ 
Whereas every entangled state is squeezed, Proposition 1 
implies that conversely any squeezed state can be entangled by 
using passive optical elements supplemented by a single addi- 
tional vacuum mode (empty port of a beam splitter), because 
the joint covariance matrix T © 1 2 [10] then satisfies inequal- 
ity (4). The optimal entangling procedure consists then of two 
steps: (i) One first applies a passive transformation S such that 
the smallest eigenvalue Ai of T is also the smallest eigenvalue 
of the 2x2 principal submatrix of S T TS corresponding to 
the first mode, (ii) One then applies the optimal entangling 
procedure on this mode and the vacuum mode, which will be 
derived in Proposition 3. 

The proof of proposition 1 leads to a lower bound for the 
attainable entanglement measured in terms of the logarithmic 
negativity [4]. The latter is so far the only calculable entan- 
glement measure for mixed Gaussian states. For an n-mode 
Gaussian state p it is given by 



Em = -^min(0,log 2 (si)). 



(15) 



where the Sj, i = 1, . . . , n are the symplectic eigenvalues of 
the partially transposed covariance matrix. Since v = \\ A 2 is 
the square of the smallest symplectic eigenvalue, we obtain 



E M > max(0,-log 2 (A 1 A 2 )/2) 



(16) 



for the attainable entanglement, with equality if there is only 
one Si smaller than one. A particularly transparent situation 
is now the case where we consider only the entanglement 
present in an arbitrary two-mode subsystem obtained when 
tracing out the other modes at the end. 

Proposition 2 Let V i — ► V = K T TK be a passive trans- 
formation acting on a Gaussian state of n > 2 modes with 
covariance matrix T. The maximum attainable amount of en- 
tanglement obtained for an arbitrary two-mode subsystem of 
r' is then given by 

E Ar = max[0,-log 2 (AiA 2 )/2], (17) 

where Ai, A 2 are the two smallest eigenvalues ofY. 

Proof: First note that for the case of a two-mode state only 
one of the two symplectic eigenvalues Si,s 2 of the partially 
transposed covariance matrix T'^ can be smaller than one, 

since (sis 2 ) 2 = dctf^ = detT'^ > 1 [13]. Following 
the same argument as in Proposition 1, there exists always a 
passive transformation S such that a two-mode principal sub- 
matrix of the covariance matrix S T TS has the two smallest 
eigenvalues Ai and A 2 , which leads to equality in Eq. (16). ■ 
A special instance of Proposition 2 is the case where the 
input already is a two-mode system, i.e. n = 2. For this case 
we will now explicitly construct the optimal entangling pro- 
cedure. We will show that it is always sufficient to perform 
a single phase rotation in one of the two modes, for example 
in A, succeeded by a beam splitter operation on both modes. 
Again, it is most convenient to employ the complex version of 
the problem. In their complex forms, a beam splitter 5(7) and 
a phase shift L(a) in system A are represented by the matrices 



S( 7 ) 



cos 7 
sin 7 



- sin 7 

cos 7 



L(a) = 







(18) 



where 7 g [0, 2tt) determines the transmission coefficient of 
the beam splitter, and a G [0, n) is the phase difference of the 
incoming and outgoing fields. Without loss of generality the 
beam splitter itself is assumed to induce no phase difference. 

Proposition 3 Let p be a Gaussian 1 x 1-mode state with co- 
variance matrix T. Let |Ai) and |A 2 ) be the eigenvectors of 
the two smallest eigenvalues ofY, with complex versions \^\) 
and |5' 2 ). The optimal entangling operation using only pas- 
sive optical elements is given by a phase rotation L(a) on 
mode A, followed by a beam-splitter 5(7/2), such that 7 and 
a are the solutions of 

cos(7)=Im[(* 2 |a z |*i)], (19) 

sin(a) sin( 7 ) = Iin[<* 2 |(7 y |*i)], (20) 

cos(a)sin(7)=Im[<*2kx|*i>], (21) 

where <j x , a y , and a z are the Pauli spin matrices. 
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Proof: In order to find the optimal entangling procedure one 
has to identify a unitary V such that 

Im [(* 2 |V£W|*i)] = 1, (22) 

and decompose it into a beam splitter and a phase shift. The 
most general form for VEV^ =: F is given by 

Fh a)=( C0S(7) Sin(7) ^ (23) 

* l7 ' } \ e M sin( 7 ) - cos( 7 ) ) ' 

which corresponds to V = L(a)B(j/2). Inserting the 
decomposition F( 7 , a) = cos(j)<t z + cos(a) sm(j)a x + 
sin(a) sin( 7 )<7 y into Eq. (23) one verifies that values a, 7 that 
satisfy Eqs. (19-21) provide a solution of Eq. (23). Moreover, 
the set of equations (19-21) always has a solution, since the 
vector of the imaginary parts in Eqs. (19-21) can be shown to 

be a unit vector if Re[(* 2 l^i)] = 0. ■ 
We will in the following apply this result to some special 
cases. The covariance matrix T of the initial state of the lxl- 
mode system will be written in the block form 

r-(£2), (24) 

where A and B are the reduced covariance matrices corre- 
sponding to mode A and B respectively. Depending on the 
form of the 2x2- matrices A, B, and C several optimal en- 
tangling protocols can be identified: 

(i) A product of arbitrary single mode Gaussian states: If 
C = 0, then a 50 : 50 beam-splitter is required in the opti- 
mal entangling procedure. The phase transformation that is 
needed will in general depend on the actual form of A and B. 
In particular: 

(ia) A product of two identical single mode states: In this 
case A = B and C = 0, and one finds that a = 7 = tt/2. The 
optimal entangling operation is thus a 50 : 50 beam-splitter, 
which follows a a — tt/2 phase transformation, as expected. 
This is the optimal procedure for uncorrected identical Gaus- 
sian input states used in several experiments [2]. 

(ib) A product of a Gaussian single mode state and a co- 
herent or thermal state: In this case where B — 61, b > 1, 
and C = the optimal entangling operation is again the ap- 
plication of a 50 : 50 beam-splitter. No phase transformation 
is required. 

(ii) States with covariance matrix in Simon normal form 
[11]: If A = at, B = bt, C = diag(c, d), then one eigenvec- 
tor is real and the other is imaginary. Hence a <G [0, 7r], 
whereas the optimal beam splitter is in general not balanced. 

(iia) Symmetric states: These are states with identical ther- 
mal reductions, meaning that A = B = at, a > 1. 
These states are already optimally entangled, since Ej^ (p) = 
max [0, — log(AiA 2 )/2], and the optimal entangling proce- 
dure is thus the identity operation. 

(iib) Special cases of symmetric states are two-mode 
squeezed pure Gaussian states with covariance matrix in Si- 
mon normal form, where in addition, C takes the form C — 
diag(c, -c) with c = (1 - a 2 ) 1 / 2 . 



In this letter we have investigated the entangling capabili- 
ties of passive optical elements in a general setting. We have 
presented a necessary and sufficient criterion for the possibil- 
ity of creating distillable entanglement in a multi-mode sys- 
tem that has been prepared in a Gaussian state. The findings 
reveal in fact a surprisingly simple close relationship between 
squeezing and attainable entanglement. We have moreover 
quantified the maximal degree of entanglement that can be 
achieved in a two-mode subsystem, and we have identified the 
optimal entangling procedure for the case of two input modes. 
In view of recently proposed applications of quantum infor- 
mation science, we hope that the presented results as well as 
the employed techniques may prove useful tools in the study 
of feasible sources of continuous-variable entanglement. 
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